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I?  ABSTRACT 
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This  paper  shows  that  there  is  a  finite  set  of  equivalence 
classes  for  constraints  in  the  0-1  programming  problem.  These 
equivalence  classes  have  the  property  that  exactly  the  same  set 
of  solutions  are  feasible  for  all  constraints  in  the  equivalence 
class.  It  is  shown  that  these  classes  are  determined  by  the  rela¬ 
tionship  of  the  constraint  to  the  dual  of  the  hypercube.  A  func¬ 
tion  that  indicates  feasibility  of  0-1  points  is  defined  and  is 
shown  to  be  monotone  over  a  vector  partial  ordering  associated 
with  the  constraint  and  paired  vertices  of  the  hypercube.  This 
allows  for  determination  of  equivalence  classes  based  on  identifying 
where  the  indicator  function  changes  value.  A  search  algorithm 
is  presented  for  classifying  the  constraints  and  a  method  is  pre¬ 
sented  for  determining  a  "best”  constraint  from  the  class. 


Constraint  Classification 
Equivalence  Classes 
Unit  Hypercube 
0-1  Programming 
Vector  Partial  Ordering 


Constraint  Classification  On  the  Unit  Hypercube. 

Abstract 

This  paper  shows  that  there  is  a  finite  set  of  equivalence  classes 
for  constraints  in  the  0-1  programming  problem.  These  equivalence 
classes  have  the  property  that  exactly  the  same  set  of  solutions  are 
feasible  for  all  constraints  in  the  equivalence  class.  It  is  shown  that 
these  classes  are  determined  by  the  relationship  of  the  constraint  to  the 

dual  of  the  hypercube.  A  function  that  indicates  feasibility  of  0-1  points 
is  defined  and  is  shown  to  be  monotone  over  a  vector  partial  ordering  associated 
with  the  constraint  and  paired  vertices  of  the  hypercube.  This  allows  for  deter¬ 
mination  of  equivalence  classes  based  on  identifying  where  the  indicator  function 
changes  value.  A  search  algorithm  is  presented  for  classifying  the  con¬ 
straints  and  a  method  is  presented  for  determining  a  "best"  constraint 


I.  Preliminaries 


Let  us  consider  two  constraints  Sa.x!  ^  b  and  Sa'x1  2b' 

1  i  ii 

where  x!  =  0, 1  Vi,  It  is  obvious  that  while  a  =  (a, ....  a  )  and 
l  In 

a'  =  (a'j, . . .  aM  may  be  outwardly  quite  different.they  may  be  redundant 
in  the  sense  that  the  feasible  sets  for  both  constraints  may  be  the  same. 

This  then  raises  the  questions  of  how  do  we  recognize  these  constraints, 
how  many  non- redundant  constraints  exist  and  how  does  one  construct  or 

choose  "best"  constraints  from  a  redundant  set.  Bales  and  Jeroslow  [4)  have 
approached  the  first  question  by  the  use  of  canonical  constraints*  The  ques¬ 
tion  of  choosing  "best"  constraints  has  been  discussed  Kianfar  [7]  and  to  some 
degree  by  Bowman  and  Nemhauser  [6J.  This  paper  addresses  itself  to  all  these 
questions  and  provides  a  means  of  classifying  constraints  into  equivalence 
classes  based  on  their  feasible  sets  and  then  choosing  a  "best"  constraint 
from  these  sets.  In  order  to  facilitate  discussion  we  shall  work  with 
variable?,  x,  that  have  been  translated  by  l/2e,  e  =  (l,...,l)'  from  the 
x1  space  i.e.  x*  =  x+I/2e*  This  translates  0,  1  values  of  x*  into 
-1/2,  +  1/2  values  of  x.  Consequently,  the  origin  in  the  x  space  is  the 
center  of  the  hypercube  whose  vertices  are  the  0,  1  solutions  desired.  We 
also  note  that  any  constraint  Sa^x^  s  bn-l/2Sa^.  In  succeeding  discussions 
we  will  be  working  only  in  the  x  space  and  thus  right  hand  side  values  of 
the  constraints  will  be  for  those  after  the  translation  into  the  x  space,  i.  e. , 
b  =  b"-  1/23^. 
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Let  K  denote  the  unit  hypercube,  B  the  minimum  ball  containing 
K,  and  C  the  minimum  octahedron  containing  K,  i.  e. , 

K  =  [x|-l/2e£x£+l/2  e]. 


For  a  detailed  discussion  of  these  sets  in  relationship  to  integer  pro¬ 
gramming,  see  f  1,2,3]. 

If  we  denote  by  bd(C)  the  boundary  of  the  set  C,  the  vertices 
of  the  hypercube  can  be  described  as  1)  bd(K)Hbd(B),  2)  bd(K)Obd(C) 
or  3)  bd(B)nbd(C). 

We  note  that  the  octahedron  H  is  the  intersection  of  2n  half¬ 
spaces 

E6.x.  sn/2  6.  =  +  1  (1) 

ii  i  — 

and  that  these  defining  hyperplanes  are  tangent  to  B  at  the  vertices  of 
K. 

Now  consider  a  hyperplane  ax  =  b  intersecting  B.  If  we  now 
construct  all  tangent  planes  to  B  at  the  intersection  Ea.x.  =  b  we 
a  displaced  cone  with  vertex  v  at  the  intersection  of  the  normal  to 
Ea.x.  =  b  and  the  tangent  hyperplanes  i.e.,  v;\a  where  X  . 
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In  projective  geometry,  v  is  the  pole  of  the  polar  set  ax  ib  with  respect 
to  the  conic  B.  We  shall  therefore  refer  to  v  as  a  pole.  When  referring 
to  the  pole  of  a  particular  constraint  ax  ab  we  shall  use  the  notation  v 
is  the  pole  of  (a,b).  For  other  uses  of  poles  and  polar  sets  in  integer 
programming  see  [2].  It  is  the  pole  v  and  its  relationship  to  the  ex¬ 
tended-facets  of  C  that  will  determine  the  equivalence  class  of  a  constraint. 

We  note  that  if  b  =  0  tlen  A  =  “  and  the  tangent  planes  form  a 
hypercylinder  (i.e.,  have  .  pole  at  infinity).  This  occurs  only  when 

ax  =  b  bisects  B,  .that. is, passes  through  the  origin,  x=0.  When  b/0 
the  pole  v  may  lie  in  either  the  feasible  closed  halfspace  ax  s  b  called  H+ 
or  the  infeasible  open  halfspace  ax  <  b  called  H~.  One  may  consider 
the  problem  of  choosing  a  "best"  constraint  as  "pushing"  the  constraint 

away  from  v  when  v  €  H“  or  "pulling"  the  constraint  towards  v  when 
v  €  H+. 

2.  Equivalence  Classes 

Let  us  now  consider  the  2n  hyperplanes  (1)  that  are  the  extended 
facets  of  C.  We  note  that  these  facets  occur  in  parallel  pairs,  i.e. , 

L6.x.  £  n/2  (2a) 

l  i  ' 

and  E(-$.)x.  5  n/2  (2b) 

ii 

I 

are  parallel.  There  are  obviously  2  of  these  parallel  pairs.  Every 
set  of  the  parallel  pairs  partition  Rn  into  3  sectors  i.e.,  feasible  to 
(2a)  but  not  (2b),  feasible  to  (2a)  and  (2b)  and  feasible  to  (2b)  but  not  (2a). 
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The  relationship  of  the  pole  of  a  constraint  to  the  parallel  pairs  (2a)  and 
(2b)  will  be  fundamental  to  the  establishment  of  equivalence  classes  on 
the  constraints.  In  order  to  discuss  pairs  we  shall  need  a  unique  designator 
for  each  pair.  Let  6  be  the  unique  designator  of  the  parallel  pair  6  and 
-6.  For  uniqueness  6  =  (6.,  6  ,...,6  )  has  the  property  E  6.  £  0 

*  m  IX  i 

and  if  E6.  =  0  then  6,  =  -1.  We  now  wish  to  establish  an  indicator 

l  1 

that  relates  the  feasibility  of  the  pole  v  of  (a,b)  to  the  constraints  (2a) 
and  (2b).  This  function  will  be  defined  over  the  2n~*  parallel  pairs 
denoted  by  6. 

Let  r(v,  6)  be  a  function  on  parallel  nairs  and  v  the  pole  of 
(a,  b)  that  exhibits  the  following  properties: 


F(v,  6) 


f  1 

if 

E6.v.  >  ~ 

xi  2 

1 

if 

S6.v.  =  ~ 
ix  2 

and  v  €  H+ 

0 

if 

E6.v.  =  j 

ii  2 

and  vCH 

0 

\ 

if 

<  E 6.v. 

2  ii 

2 

1 

\  0 

1 

if 

S5.V.  =  -f 

ii  2 

and  vc H 

i 

-i 

if 

S6.v.  =-f 

ii  2 

and  v  €  H+ 

-i 

V 

if 

E6.v.  <-§ 

ii  2 

T<v,  5)  describes  in  which  of  the  three  partitions  the  point  v  lies.  We 
thus  have  the  following  Lemma. 

The  extended  facets  of  the  octahedron  C  partition  Rn  into 


Lemma  1: 


■5 


Proof:  There  are  2  parallel  pairs  6  and  any  point  v  is  mapped  into 


one  of  the  three  possible  values  under  F(v,  6)  depending  on  where  v  lies 


in  relation  to  the  partition  of  Rn  by  the  pair  6.  There  are  2n~*  such 


mappings  for  a  given  v . 


We  now  investigate  the  relationship  of  the  pole  v 


of  the  hyperplane  ax  =  b,  to  feasible  vertices  of  K  in  the  halfspace  ax  ^b. 


Since 


the  pole  v  is  said  to  have  X  =  “  for  b  =  0,  we  extend  our  definition  of 


r  to  include  these  values. 


1  if  a6<  0 


J7(®,  6)  =  0  if  afi  =  0 


-1  if  a&>  0  . 


Recall  that  1/26  and  -1/26  are  vertices  of  K,  i.  e. ,  solutions  to  our 


problem.  The  following  Lemma  relates  the  feasibility  of  these  vertices 


to  the  function  T(v,  6 ). 


Lemma  2:  Let  ax  ^b  have  a  pole  v  then 


if  b  >  0 


a)  1/2  6  is  feasible  if  and  only  if  T(v  ,  6)  =  +1 


b)  -1/2  6  is  feasible  if  and  only  if  T(v  ,  6)  =  -1 


if  b  <  0 


c)  1/2  Sis  irfeasible  if  and  only  if  r(v,  6)  =  1 


d)  -1/2  &s  infeasible  if  and  only  if  T(v,  6)  =  -1 


and  if  b  =  0 


e)  1/2  &s  infeasible  if  and  only  if  r(°°,  6)  =  1  i.e.,  a6<0 


f)  -1/2  Sis  infeasible  if  and  only  if  rf”,  6)  =  -1  i.  e. ,  a6>0 


l.  e 


Proof: 

a)  Assume  1/2  6  is  feasible,  then  1/2  6a  ^  b  or  multiplying  by  “£■ 

we  have  ~”~a6  2  ~  or  v6  2  Since  b>  0  then  v€H*  and  thus 
4d  i  z 

T(v,  5)  =  1.  Now  let  T(v,  6)  =  1  then  v6a~  or  ~  a  6  2  ^  or 

2  4b  c 

1/2  a  6  2  b.  The  proof  of(b)  is  similar. 

To  prove  (c)  let  1/2  6  be  infeasible  then  a(i/2  6)  <  b  or  multiplying 

by  ~  <  0  we  get  (~)  a  5  >  ~  and  thus  T{v,  6)  =  1,  Now  assume 
2b  4b  L 

that  r(v,  6)  =  1.  Since  b  <  0  then  vCH  and  we  have  (*~)afi  >  ~ 

4b  2 

and  consequently  1/2  a6<  b.  Again  the  proof  of  (d)  is  similar  to  that  of 
(c).  To  prove  (e)  assume  1/2  6  is  infeasible  then  a(l/2)6  K  0  or. 

a  5  <  0.  Likewise  if  a  6  <  0  then  1/2  a  6  <  0.  Again  the  proof 
of  (f)  is  similar  to  that  of  (e). 

We  have  thus  shown  that  we  can  use  the  location  of  the  pole, 
of  (a,  b)to  yield  information  about  feasible  and  infeasible  vertices  of  K 
with  respect  to  ax  2  b. 

Definition.  We  say  that  two  constraints  (a,  b)  and  (a',  b’)  are  in 
the  same  equivalence  class  E,  if  the  set  of  feasible  vertices  of  K  are 
the  same  for  both  (a,  b)  and  (a1,  bf). 

We  denote  (a,  b)  in  equivalence  class  E  by  (a,  b)  €E. 

^n- 1  2n" 1 

Theorem  1.  There  are  (2*  3  -  2  )  equivalence  classes  for 

constraints  in  an  n-dimensional  hypercube. 


Proof:  From  Lemma  2  we  can  classify  opposite  points  1/2  0  and 

-1/26  by  the  function  T.  Thus 
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a)  1/2  6  and  -1/26  are  both  feasible  if  and  only  if  b  £  0 

and  T(v,  5)  =  0 

b)  1/2  f  is  feasible  and  -1/26  are  infeasible  if  and  only  if 

1)  b  >  0  and  T(v,  6)  =  1 

or  2)  b  s  0  and  T(v,  6)  =  -1 

c)  1/2  6  is  infeasible  and  -1/26  is  feasible  if  and  only  if 

1)  b  >  0  and  T(v,  6)  =  -1 

or  2)  b  s  0  and  T(v,  6)  =  1 

d)  1/2  6  and  -1/26  are  both  infeasible  if  and  only  if 

b  >  0  and  r(v,  6)  =  0 


We  first  note  that  there  are  3  values  of  T(v,  6)  for  each  of  the 
1  — 

2  parallel  pairs  6  and  that  the  feasible  sets  differ  with  the  sign  of  b  (that  is 

b  £  0  or  b  >0).  Thus  an  upper  bound  on  the  number  of  equivalence 
-n-  1 

classes  is  2(3^  ).  However  (b)  and  (c)  above  indicate  that  different 

2n“  * 

elements  of  these  2(3£  )  potential  classes  are  the  same.  In  particular 

if  r(v,  5)  =  +  l  then  changes  of  the  sign  of  T  and  b  produce  the 
same  equivalence  class.  Thus,  since  there  are  2  ways  of  assigning 
f{v,  6)  =  +1  for  all  6  and  2  ways  of  assigning  the  sign  of  b,  there 

•tial  classes  of  which  1/2  are 
the  same  i.e.,  the  change  of  sign  of  T  and  b.  Thus  from  the  total 
set  of  potential  classes  we  must  subtract  those  that  are  redundant.  This 
yields 


.  .  ...  ,  jz"-1  +  1) 

is  a  totality  of  2  of  thes 


_n- 1 


.n-1 


2(3 


)  -  2 


equivalence  classes. 
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I 

& 

! 

It  is  obvious  that  these  equivalence  classes  can  be  determined  by  } 

! 

just  observing  the  values  of  T  and  b  as  outlined  in  (a)  through  (d) 

I 

of  the  above  proof. 

Example  li  Let  us  consider  the  possible  equivalence  classes  in  a  two 
dimensional  space.  There  are  four  vertices  of  K  which  we  denote 
U,  X.  Y,  Z.  We  let  U  =  (-1/2,  -1/2),  X  =  (1/2,  1/2),  Y  =  (-1/2,  1/2) 
and  Z  =  (1/2,  -1/2)  thus  these  correspond  to  the  following  0-1  points 
before  translation  of  axis*  U'  =  (0,  0),  X'=  (1,  1),  Y*=  (0,  1)  and  Z‘=  (1,  0). 

Let  7) 1  =  (1, 1)  and  5  2  =  (-1, 1)  then  U,  X  represent  the  vertices 
classifiable  by  r(v,6*),  and  Y  and  Z  the  vertices  classifiable  by 

62).  The  number  of  possible  equivalence  classes  is  determined  from 

1  2  2 
Theorem  1.  There  are  2n“  =  2  pairs  of  vertices  and  thus  2»  3  -2  =  14 

possible  equivalence  classes.  These  are  as  follows 


b 


1) 

*0 

2) 

SO 

3) 

SO 

4) 

so 

5) 

SO 

/  t 

rso 

7) 

\>l 

f  so 

8)  . 

\>o 

i 

/so 

9) 

{>0 

10) 

>0 

ID 

>0 

12) 

>0 

13) 

>0 

14) 

>0 

r(v,  61) 


r(v,  a2) 


feasible  vertices 


0 

0 

0 

-1 

1 

-1 

1 

-1 

1 

1 

-1 

1 

-1 

0 

0 

0 

-1 


U,  X,  Y,  Z 
U,  X,  Y 
U,  X,  Z 
x,  Z,  Y 
U,  Z,  Y 

<x,  v> 
<X.  z> 
{u,  y> 
{u,  z> 

none 

Z 


1 


Y 


0 


U 


1 


0 


X 


I 


1 


I 
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One  should  note  that  the  number  of  equivalence  classes  grows  quite 
rapidly  as  shown  below  !  1 


n 


Equivalence  classes' 


2 

t 

14 


3  .  . 
'146  ‘ 


4 

12866 


86,027, 906 


While  we  have  been  able  to  classify  constraints  into  equivalence  , 

;  I  1 

classes  by  Theorem  1,  for  any  given  constraint  we  must  calculate  2 

i 

functions  F(x,  6)  *  to  determine  this  classification.  Although  this  is  not 

l  % 

equivalent  to  enumerating  all  of  the  vertices  of  the  hypercube,  it  is 

1  \  ) 

equivalent  to  enumerating  half  of  the  vertices.  In  the  next  section  we  - 

1  i  1 

investigate  a  partial  ordering  that  exists  on  the  extended  ‘facets  of  C 

I  ;  ' 

and  from  this  derive  an  efficient  search  algorithm  for  determining  the 

l 

.  1  , 

appropriate  equivalence  class.  In  addition  th,is  partial  ordering 'suggests 

•  -  i  s  . , 

a  simple  method  of  deriving  other  constraints  in  the  equivalence  clads 


that  are  "better". 


i  * 


i 


1 
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3.  Calculations  of  T(* ,  •) 

If  we  observe  that  the  hypercube  may  be  rotated  and  the  axes 
relabeled  we  can  assume  without  loss  of  generality  that  ax  2b  has  the 
property  0  s  aj  ^  a^  s  s  a^.  Since  the  6  for  paired  hyperplanes  has 
the  property  26.  2  0  we  conclude  that  at  most  [rr  J  of  the  6  can  be  negative 
([x]  is  the  greatest  integer  s  x).  We  can  thus  describe  6  in  terms  of 
the  negative  elements  of  6. 

Let  k  58  [n/2],  p  =  be  a  Ic-dimensional  vector,  and  k*(t) 

be  the  number  of  negative  elements  of  3,  With  every  6  we  shall  associate 
a  unique  vector  p(6)  that  has  the  following  properties: 

a)  0  <  p^  <  Pj,  <  .  • .  <  p^ 

b)  if  pt  4  0,  then  Pi  <  P1+1  (3  l) 

O  Pt  -  P2  -  ...  =  Pk.k.(5)  -  0. 

d)  6^  *=  -1  if  and  only  if  p^  =  i  for  some  j. 

Property  (d)  establishes  the  relationship  for  identifying  6  given  p(5), 
i.e.,  p(J)  is  a  list  of  the  indices  of  the  negative  elements  of  6.  The 
upiqueness  of  p(6)  is  established  by  properties  of  (a) ,  (b)  and  (c)  •  For  - 

notations 1  convenience,  we  shall  use  the  notation  p  for  p(6) . 

Example  2:  For  n  =  5  k  =  2  some  examples  of  associated  p  are 


i 


6  =  (1,  1,  1,  l,  1) 

6  =  (-1,  1,  1,  1,  1) 

6  =  (1,  -1,  1,  -1,  1) 


P  =  (0,  0) 
P  =  (0,  1) 
P  =  (2,  4) 
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The  p  values  introduce  a  vector  partial  ordering  where  p  >  p*  it  and 
only  if  p^  >  for  all  i.  If  strict  inequality  must  hold  in  at  least  one 
component  we  write  p  >  p*.  This  partial  ordering  f,a  important  with  res¬ 
pect  to  the  indicator  function  r(*,*)>  thus,  in  order  to  maintain  notational 
consistency,  we  define  a  function  H(p;v)  as  follows: 

m  *  '*  ■**'  ’  m 

Definition:  H(p;v)  *  T(v,6)  where  p  is  the  vector  associated  with  6. 
Definition:  Two  vectors  p  <  p*  are  adjacent  if  ji^  =  p^|  for  j  4  i  and 
p*  *  p±  +1* 

We  can  now  establish  an  important  roonotonicity  relation  on  the  function  H. 
Theorem  2:  If  p  <  p'  then  H(p;v)  >  H(p' ;v)  for  b  >  0  and  H(p;v)  <  H(p*;v)  for 

. ■■■  hi  ■— — i — w  i  — — — —— ■» 

b  <  0. 

Proof:  Consider  p  and  p*  such  that  p  <  p*  and  p*  is  adjacent  to  p.  Let 
p*  *  p  +  1.  The  condition  p  <  p*  and  property  (a)  of  (3.1) 

Pj  4  pi  +  1  for  any  j  since  otherwise,* he  number  of  non-zero  values  p* 
would  be  one  less  than  that  of  p  by  property  (b)  of  (3.1).  We  also  note 


and 
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X£  6  (p)a  =  X  £  8  (p)a  +X(a  -  a  ) 

t=l  ‘  *  «q,q+l  ‘  ‘  q+1  15 


Since  a  .  -  a  ^  0  by  assumption  then  H(p;  v)  «s  H(p*;  v)  if  X  >  0 
q+l  q 

i.e.  b>0  and  H(p;  v)  £  H(p*,  v)  if  X  <  0  i.  e.  b^O. 

* 

We  now  need  to  show  that  for  any  p  ^  p',  p*  can  be  reached  by  a 
path  through  adjacent  vectors.  Let  =  q  and  p^  =  q  +. r.  Then  let 


p°  =  p  and  pO  5  p1  s , . ,  s  pT  where  p*  =  p*”*  +  1  for  i=l,...,r 

K  xv 

and  p.  =  p.~  for  j  ^  k  .  This  establishes  montonicity  of  H  on  p°  through 
J  3 

pr  by  the  above  proof  on  adjacent  vectors  p1  and  p  .  Also,  p1  <  p'-.  We  oow 

take  index  k-.l  and  perforin  the  sane  operation  and  continue  with  index  k-2,  k-3, 

...,  1  and.  thus  have  p  <  p°  <  p*  <  ...  <  pm  *  p*  where  p**^  and  p*  differ  by 

only  one  element.  Thus,  the  monotonicity  property  is  preserved  and  the  theorem 
is  proved. 

Theorem  2  provides  the  basis  for  calculating  the  T(v,  6)  by  a 
search  algorithm.  Note  that  by  the  theorem  one  has  only  to  calculate  the 

vectors  p  where  the  function  H  changes  value.  Thus,  if  H(p)  *  1  and 

H(p*)  =  0  with  p  and  p*  adjacent  and'  b  >  0,  we  have  H(p*)  =  1  for 

p'  <  p  and  H(p')  =  0  or-1  for  p*  >  p*,  Thus,  if  we  describe 

the  set  of  p  where  H  changes  values  we  completely  describe  the  function  T 

and  thus  the  equivalence  class  in  which  (a,b)  exists.  Towards  this  end, 

we  define  the  following  sets  that  characterize  these  points  and  thus, 

the  equivalence  class. 
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If  b  *  0  then 

P+  =  {  p|H(p;v)  =  1  and  H(p*;v)  *  0  or  -1,  p*  <  p} 

P0+=  {p  |H(p;v)  =  0  or  -1  and  H(p*;v)  »  +1  for  p*  >  p} 

P^  =  { p  |  H(p;v)  =  0  or  +1  hnd  H(p'jv)  =  -1  for  p'  <  p} 
and  P*  =  {p|H{p;v)  =  -1  and  H(p*;v)  =  0  or  1,  p*  >  p} 

and  if  b  >  0  then 

P+  =  £  p  J  H(p;v)  =  1  and  H(p‘;v)  =  0  or  1,  p*  >  p} 

P^+=  j-p|H(p;v)  =0  or  -1  and  H(p*»v)  -  +1  p*  <  p} 

P^~  =  {p  j  H(p;v)  =  0  or  +1  and  H(p*;v)  =  4  p'  >  p} 
and  P~  -  £p|H(p;v)  =  -1  and  H(p';v)  =  0  or  1  p*  <  p} 

From  Lemma  2,  the  sets  P+  and  p“  provide  characterization  of  the  feas¬ 
ible  points  as  follows: 

a)  1/2?  and  -1/26  are  both  feasible  if  and  only  if  b  <  0  and 
for  some  p*  c  P  and  some  p”  6  P  ,  p*  <  p(5)  5  PM*  (3»2a) 

b)  1/26  is  feasible  and  -1/26  is  infeasible  if  and  only  if  either 

1)  b  >  0  and  for  some  p'  €  P+  p(6)  <  P*  (3.2b) 

or 

2)  b  <  0  and  for  some  p'  €  P  p(6)  <  P* . 

c)  1/26  is  infeasible  and  -1/26  is  feasible  if  and  only  if  either 

1)  b  >  0  and  for  some  p'  €  P  p(6)  >  p*  (3.2c) 

or 
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d)  1/25  and  -1/26  are  both  infeasible  if  and  only  if  b  >  0  and  for  aome 
p*  6  P+  and  p"  6  P“ 

P*  <  P(6)  <  p".  <3. 2d) 

The  sets  P^  and  P®  are  used  in  the  algorithm  of  the  next  section  to  de- 
— 

termine  P  and  P  and  will  be  used  to  generate  linear  constraints  des¬ 
cribing  the  appropriate  equivalence  class  in  Section  5*.  ..  . 


Example  3:  Consider  the  constraint 

2x,  +  4x_  +  5x,  +  7x  +  8x_  £  -4 
1  2  3  4  5 

There  are  16  possible  vectors  p.  The  pole  is  v  =  -5/16(2,  4,  5, -7,  8) 
and  the  values  of  H  are  as  follows: 


p 

6(p)v 

H(p,  v) 

P 

6(p)v 

00 

65 
"  8 

-1 

14 

20 
••  IT  ir 

8 

01 

55 
"  8 

-1 

15 

15 
"  8 

02 

45 
~  8 

-1 

23 

20 
"  8 

03 

«|o 

-1 

24 

10 
’  8 

04 

30 
'  8 

-1 

25 

5 

8 

05 

25 
*  8 

-1 

34 

5 

“  8 

12 

35 

8 

-1 

35 

0 

13 

30 
‘  8 

-1 

45 

10 

8 

H(p,  v) 
0 
0 
0 
0 
0 
0 
0 
0 


Now  we  observe  that  H(0,5;v)  =  -1  implies  H(0,t;v)=-1  for 
t  =0,1...,  5  and  H(l,  5;v)  =  0  implies  H(j,  5;  v)  =  0,  or  1  for 
j  =  1,  2,  3,  4,  and  since  H(4,5;  v)  =  0  then  these  are  all  zero.  Since 
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all  of  the  values  by  H(0,5)  -  -1,  H(l»3)  =  -1,  H(lj4)  *  0  and  H(2,3)  **  0« 
Alternatively,  we  have  P*  *  {(0,5),  (1,3)},  P°  =  £(1,4) , (2,3)} ,  ** 


In  the  next  section  we  describe  an  algorithm  for  finding  the  seta 


P. 
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4.  An  Algorithm 


In  this  section  we  present  a  method  for  determining  the  sets  P  , 

Of-  0- 

P  ,  P  and  P  .  The  method  determines  elements  in  one  set  and  then 
insures  that  all  adjacent  elements  are  dominated  (in  a  partial  ordered 
sense)  by  the  elements  of  one  of  the  sets.  In  the  algorithm,  we  shall 
be  interested  in  those  adjacent  vectors  that  are  either  greater  than  or 

less  than  a  particular  vector.  Towards  this  end  we  use  the  following 

notation.  The  vectors  adjacent  to  p*  will  be  denoted  as  p* ( j)  j  *  l,...,n^ 

and  Pg( j)  j  =  l,...,n£  where  p*  ( j)  <  p'  for  all  j  and  p^(j)  >  P*  for 
all  j.  Thus,  the  p[(j)  are  the  adjacent  vectors  to  p*  that  are  less 
than  p*  in  partial  ordering  and  p^( j)  are  the  adjacent  vectors  that  are 
greater  than  p' . 

In  addition,  since  the  monotonicity  of  H(p}v)  is  reversed  for  rever¬ 
sal  of  the  sign  of  b  the  algorithm  is  expressed  for  b  >  0  with  changes 
for  b  <  0  designated  in  brackets.  Finally,  since  we  are  constructing 
sets,  we  use  +  and  -  to  denote  the  addition  and  subtraction  of  elements 
from  sets. 

Step  0:  P+  =  P0*-  =  P~  =  0.  T+  =  T0*"  =  T°"  =  T*  *  0. 

Step  1:  Initial  Search:  Search  an  ordered  chain  from  p  »  (000,..., 0)  to 

p  =  (p^,...,n-l,n),  where  p^  =  n-k-1  if  n  is  odd  and  pj»  1  if  n  is  even,  let 
p*  be  the  maximum  element  such  that  H(p*,v)  =  1  [H(p*,v)  =-  -lj  and  p**  be 
the  minimum  element  such  that  H(p**,v)  =  -1  [H(p**,v)  =,  +11.  ■  Let  T  -  {p,*} 
IT~  -  {p*}]  and  T*  =  {p**}  [T+  =  {p**}3. 

Remark:  It  is  possible  that  p*  and  p**  may  not  exist*  In  this  case,  the 
appropriate  set  is  left  null.  While  any  ordered  chain  may  be  used,  a 


PilL 


convenient  one  is  (00,...0,j)  j  =  0,1,. ..,n  ;  (000,. . . ,0, j,n)  j  =  l,...,n-l  J 
(000,.. .  ,0,j,n-l,n)  j  =  l,...,n**2  ,  etc. 

Step  2:  If  T+  =  0  go  to  Step  3;  otherwise,  select  any  p1  €  T+  and  T+  «* 

T+  -  £p'}. 

4- 

2a:  If  p'  <  p  [p^>  pj  for  any  p  €  P  ,  go  to  Step  2;  otherwise, 
j  =  1  go  to  2b. 


2b:  If  H(p'(j),v)  =  1  tH(p*(j>;V)  =  1],  let  p*  =  pj(j)  [p*  *  p'(j)l 
go  to  2a; 

otherwise,  j  =  j  -f  1  and  go  to  2c. 


2V2 


2c: 

If  j  >  n’  [j  >  n^J  go 

to  2d; 

2d: 

P+  -  p++  {p‘3.  ^  = 

-f 

[T0f  =  Ttx-  +  {p»  (1)]  . 

f  . . »  + 

Step  3: 

Of 

If  T  =0,  go  to  step  4. 

Otherwise,  select  any 

p'  €  T* 

3a: 

If  p'  >  p  Cp*  <  p]  for  any  p 

Otherwise,  j  =  1,  go 

to  3b, 

3b: 

If  H(p‘(j)i;v)  =  0  or 

-l[H(p' 

p{(j)lp*  =  P2(j)3  go 

to  3a;  i 

3c: 

If  j  >  tj  >  n^J  go 

to  3d; 

Otherwise,  go  to  3b. 

■ 

3d: 

P*  -  P0”  +  (p'5.  T+ 

=  T+  + 

it+  =  t+  +  {,;<!)}+  . 

•  •  +  £p. 

Of  .  ^Of 
and  T 

„0f  __ 


„0f 


*f  Of 

Remark:  Upon  entering  step  4  we  have  determined  the  sets  P  and  P 
Step  4  and  5  are  analogous  to  Step  2  and  3  and  determine  P  and  P* 
Step  4:  If  T*  =  0,  go  to  Step  5;  Otherwise,  select  any  p'  €  T*  and 
T~  =  T”  -  {p*3 . 


sagfia 
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4a: 

If  P1  >  P  fp*  <  p] 

for  any  p  €  p“,  go 

to  Step  4;  Otherwise, 

j  *  1,  go  to  4b. 

4b: 

If  H(p»(j)‘,v)  =  -1 

fH(p2{j)»v)  =  -1]  let  p'  =  p^(j) 

fp*  =  p'U>] 

go  to  4a;  Otherwise,  j  =  j  +  1  and  go 

to  4c. 

4c: 

If  j  >  n*  [j  >  n*] 

go  to  4d;  Otherwise 

,  go  to  4b. 

4d: 

p"  -  p“+  Cp'},  t°* 

A 

’  =  I0"  +  CpJ(1)}  + 

...  +  {pj(np} 

r_0— 

[T  * 

T  “+  {p^(l)}+  ... 

'+  CpJ(nJ)}].  GO  to  4. 

Step  5: 

If  T  -  0  stop;  Otherwise,  select  any 

p’  €  T°"  and 

.  0- 
let  T  * 

t°‘  -  CP'}. 

5a: 

if  p*  <  P  [p*  >  pj 

for  any  p  €  P°“,  go 

to  Step  5:  Otherwise, 

j  =  1.  Go  to  5b. 

5b: 

If  H(p*(j)yv)  =  0  or  +1  [H(p^(j);v)  = 

0  or  +1 ] ,  let 

P'  =  Pj(j) 

tp*  “  p[(j)l.  Go  to  5a;  Otherwise,  j 

=  j  +  1  and  go 

to  5c. 

5c: 

If  j  >  n’  [J  >  njl. 

go  to  2d;  Otherwise  go  to  5b. 

5d: 

P°‘  -  p°'+  {p*J.  t 

‘  »  T*+  {p!(J)}  +  , 

•  ••  +  {po(<)3 

It"  » 

T  +  {p[(l)}  +  ...  +  ,  {p^(tt*)}3.  Go  to  4. 

In  the  next  section,  we  discuss  a  possible  application  of  equivalence 


classes  to  determine  a  "best'* 


constraint  for  the  equivalence  classes. 


5.  Generating  Constraints  from  Equivalence  Classes 


Let  us  now  suppose  that  we  have  an  equivalence  class  as  represented 
by  the  sets  P*  and  P  .  A  question  then  arises  as  to  what 
constraint(s)  should  we  use  to  represent  this  equivalence  class  and  how 
can  they  be  obtained.  We  represent  our  new  constraint  as 

^1X1  +  ^2X2  +  •**  *  ^nXn  -  ^0* 

Since  H(p;v)  is  a  function  after  both  rotation  and  permutation  of  the 
axis  have  been  eliminated  we  have  the  additional  conditions 

0  <  8,  <  8.  ...  <  8  . 

In  addition,  we  also  know  that  the  class  is  determined  by  the  sets  P+ 
and  P  of  the  function  H(p;v).  We  therefore,  construct  a  set 
of  linear  constraints  using  these  nets  and  and  P^~. 

From  the  sets  P+ ,  P^,  P^  and  P  and  (3.2)  we  have  the  fc  1  lowing 
conditions: 

If  b  <  0  then 


-l/26(p)8  <  8q  - 

1/2 

for 

p  €  P* 

(5.1a) 

1/26 (P) 8  <  8q  - 

1/2 

for 

p  6  P+ 

(5.1b) 

-l/28(p)p  >  8q 

for 

p  €  P°" 

(5.1c) 

l/26(p)P  >  8q 

for 

P  €  P04- 

(5. id) 

and  if  b  >  0  then 


-l/26(p)P  >  pQ 

for 

p  €  P" 

(5.1e) 

1/26 (p) 8  >  80 

for 

p  €  P+ 

(5. If) 

-l/26(p)P  <  80  “  1/2 

for 

,  0~ 

P  €  P 

(5.1g) 

l/26(p)8  <  -  1/2 

for 

p  €  P^ 

(5.1h) 

j 
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These  conditions  come  directly  from  Lenina  2.  We  note  that  P  and 


P^"  are  those  points  that  are  adjacent  in  the  sen  se  that  if  b  <  0 ,  -l/26(p) 


is  feasible  for  p  6  P^"  while  -l/26(p)  is  infeasible  for  p  €  ?*  .  This 


fact  is  reflected  in  statement s(5.1c)  and  (5.1a).  We  use  3q~  1/2  for 


infeasibility  since  this  translates  back  to  3q  -  1  (i.e.,  integer  r.h.s.) 


indicating  infeasibility  in  the  x*  (zero-one)  space.  The  linear  con¬ 


straints  (5.1)  have  the  property  that  all  3  =  (3^,3^, ...,3^)  satisfying 


them  belong  to  the  same  equivalence  class  and  any  3  that  does  not.  satisfy 
these  constraints  does  not  belong  tA.the  same  equivalence  class 
as  the  constraint  (ax  >  b)  from  which  the  sets  P+,P  were  derived. 

The  criteria  of  choosing  a  "best"  constraint  from  the  set  of  solu¬ 
tions  to  the  linear  constraints  will  of  course  depend  on  the  definition 
of  the  term  "best".  We  propose  that  one  meaningful  definition  is  that 
3  provides  a  best  constraint  if  the  member  of  binding  (equality)  con¬ 
straints  for  the  defining  linear  inequalities  (5.1)  is  as  large  as 
possible.  From  linear  programming  theory  we  know  that  these  conditions 
corresponds-  to  generating  basic  feasible  solutions  for  the  constraints 
(5.1),  where  the  number  of  non-basic  slacks  is  as  large  as  possible. 

This  immediately  implies  that  there  may  be  several  "best"  constraints 
since  there  may  be  several  such  basic  solutions.  However,  they  all 
possess  the  property  that  a  slight  pertubation  of  a  coefficient  either 
makes  one  of  the  defining  equalities  non-binding  or  makes  3  no  longer 
satisfying  the  constraint  set.  It  should  also  be  noted  that  there  are 
nfl  variables  in  (5.1).  While  this  may  seem  disconcerting 
since  we  need  only  n  points  to  determine  a  hyperplane  in  n  space,  it  is  a 
result  of  our  definition  of  an  infesible  point  as  being  at  least  1  unit 


irm-i  fTf  n 
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length  away  from  the  constraint  in  terms  of  the  right-hand  side.  One 
may  invision  this  process  as  choosing  n  points  to  determine  the  hyper¬ 
plane  and  then  one  additional  point  so  that  a  displacement  of  the  hyper¬ 
plane  by  one  unit  (in  the  x'  space)  will  pass  through  an  infeasible  point, 
see  (5.1a,b,g  or  h). 

\ 

%*•  i'  ' 

Remark:  v  v 

We  note  that  the  constraint  (5.1)  can  be  translated  back  to  the  x*  space 
since  x^  =  -1/2  implies  =  0  and  x^  =  1/2  implies  x*  *  1#  Since  this 
reduces  the  number  of  non-zero  coefficients  in  the  constraint  set  pro¬ 
duced,  we  propose  that  this  should  be  used. 

Example  4: 

Let  us  now  consider  the  constraint  of  Example  3.  In  that  example 
we  had  =  P+  =  0,  and  p”  =  [(0,5),  (1,3)]  and  P°  =  [(1,4)  ,(2,3)  ]. 

The  constraints  from  (5.1)  are  thus, 

s5  <  s;  -  i 
B1  +  %  s  Bo  *  1 

+  S>4  ^  Bj  <5-2> 

B2+S3 

0  <  Bj  <  S2  <  S3  <  S4  <  S5 

The  first  constraint  comes  from  the  fact  that  p  =  (0,5)  implies  o  =  (1,1,1, 
1,-1)  and  since  -1/26  is  infeasible,  then  x|  =.  1,  =  x^  =  *  0 

is  infeasible.  The  other  constraints  are  derived  similarly. 

We  also  note  that  we  can  change  variables  and  eliminate  the  con¬ 
straints  o  <  <  33  <  34  <  35*  let  y^  =  9^  and  yjL  ~  -  0.^ 

for  i  >  2  then  y^  >  0  implies  0^  >  0^^  thus,  (5.2)  reduces  to  finding  a 


solution  to 


yl  +  y2  +  y3  +  y4  +  y5  -  0O  “  1 

2yl  +  y2  +  y3  ^  *'o  "  1 

2yL  +  y2  +  y3  +  y4  >  $'0 

2yi  +2y2  +  y3  2 

yt  >o 

Solving  this  set  of  constraints  we  find  that  =  1,  y^  =  1,  y^  =  0, 
y4  -  1*  y^  =  0,  (3^  =  4  is  a  basic  feasible  solution  (and  the  only 
basic  feasible  solution).  This  yields  the  constraint  x'^  +  2x^  +  2jfc^ 

3x^  +  3x^  >  4.  It  is  easy  to  check  that  this  is  from  the  same  equiva¬ 
lence  class  as  the  original  constraint. 

Remarks:  1)  The  equations  (5.1)  indicate  that  the  flexibility  of 
"moving"  constraints  and  maintaining  the  same  feasible  vertices  of  the 
hypercube  is  dependent  not  only  on  keeping  certain  vertices  feasible 
i.e.,  equations  (5.1c),  (5. Id),  (5.1e)  and  (5. If)  but  also  on  keeping 
certain  vertices  infeasible,  i.e.,  equations  (5.1a).  (5.1b),  (5.1g) 
and  (5.1h). 

2)  One  should  not  interprets  the  constraint  generation  of  this 
section  as  being  the  only  or  even  the  main  application  of  theory  of  con¬ 
straint  equivalence  classes.  In  15 J,  the  concept  of  equivalence  classes 
is  directly  applied  to  a  special  class  of  0-1  programming  problems  re¬ 
sulting  in  not  only  the  optimal  solution  but  direct  sensitivity  analysis 
on  the  changes  in  the  objective  function. 


mvctia 
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